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Isospectral Potentials from Modified Factorization 

Micheal S. BergeiB and Nail S. Ussembaye\0 
Department of Physics, Indiana University, Bloomington, IN 4^405, USA 

Factorization of quantum mechanical potentials has a long history extending back to the earliest 
days of the subject. In the present paper, the non-uniqueness of the factorization is exploited to 
derive new isospectral non-singular potentials. Many one-parameter families of potentials can be 
generated from known potentials using a factorization that involves superpotentials defined in terms 
of excited states of a potential. For these cases an operator representation is available. If ladder 
operators are known for the original potential, then a straightforward procedure exists for defining 
such operators for its isospectral partners. The generality of the method is illustrated with a number 
of examples which may have many possible applications in atomic and molecular physics. 
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PACS numbers: 03.65.Ge, 03.65.Fd, 03.65.Ca 

I. INTRODUCTION 

The factorization method due to Hull and Infeld [l[ has 
been widely exploited in quantum mechanics to deter- 
mine the spectra and wave functions of exactly solvable 
potentials. This approach has been formalized in super- 
symmetric quantum mechanics (SUSY QM) |2| which has 
been used to find many new isospectral potentials. The 
usual procedure is to find a factorization of a quantum 
mechanical Hamiltonian and the methods of SUSY QM 
then guarantee that a supersymmetric partner potential 
is isospectral to the original Hamiltonian. As verified be- 
low, this procedure yields a pair of potentials with the 
same spectra (possibly apart from the ground state) and 
related wave functions. Throughout this paper we work 
in h = 2m = 1 units. 

Let's consider a one dimensional Hamiltonian 



H 



(0) 



-d 2 x + V m (x) 



where (x) is an arbitrary non-singular potential with 
at least one bound state and zero ground state energy 
(given the Hamiltonian H = —d 2 + V(x) one simply sub- 
tracts the zero point energy to obtain H^). It is a sec- 
ond order linear operator and it can be factored into a 
product of first order linear operators as follows: 



#i 0) = (-d x + W (x))(d x + W (x)) = AlA 

once the ground state wave function ipo(x) is specified. 
The function Wq{x) = —d x \nip (x) is called superpoten- 
tial generating the potential 

V^(x) = W^x)-W^x). 

Fortunately, the factorization does not commute 
AqAq AqAq unless the superpotential is constant. In 



other words, an inverted product ^o^o ^ s a certain new 
Hamiltonian H ( ° ) = A G A\ = -d' 2 x + u| 0) (a;) where 

vi Q) (x) = W 2 (x) + W^x) 

is also free of singularities. It turns out that the eigen- 
functions and eigenvalues of these partner Hamiltonians 
are related. Indeed, we have the following first-order in- 
tertwining relations 



H^Al = Ajtff and fff A 



A oH M 



(1) 
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from which one observes that since Aoip (x) = 0, the 
spectra of H+ and are connected by E n = E n+ \ 
(n = 0,1,...) where E n and E n denote the eigenval- 
ues of the Hamiltonians and respectively with 
eigenfunctions ip n and yb n . Thus, the Hamiltonians have 
identical energy spectrum except for the ground state of 
H^°\ The wave functions satisfy ^ n {x) oc Aoip n+ i(x), 
fpn+i{x) oc A , ip n (x) and if ip n+ i(x) is normalizable, then 
ipn(%) is also normalizable and vice versa, because 

(ip n (x),ip n (x)) = (tp n+i (x),AlA il} n+1 (x)) 
= E n+ i(tjj n+ i(x) 1 ip n+1 (x)). 

Note that for singular potentials (for instance, with a 
1/x 2 singularity) some of the wave functions ip n (x) are 
not acceptable as they may not be normalizable [3| . That 
is, for singular potentials the degeneracy of energy levels 
is only partially valid or invalid at all. The upshot of all 
this is that one can generate new isospectral potentials 
from existing exactly solvable potentials. 

Luckily, the above discussed factorization is not 
unique. For example, we have 

(-d x + l){d x + I) = {-d x + ta,nh(x))(d x + tanh(x)), 

i.e. two different superpotentials can give rise to the 
same potential (in this particular example with no bound 
states). One can try to construct new isospectral poten- 
tials exploiting non-uniqueness of factorization and ob- 
tain a one-parameter family of potentials with the pa- 
rameter arising as an integration constant 0, Q • 
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Suppose the Hamiltonian can be factorized by the 
operators different than Aq and Aq, namely, 

B = 8 X + f{x) and St = -d x + f(x) 

where f(x) is temporarily undetermined function: 



ijf = BB 



t _ 



= -d 2 x + f(x) + f(x). 



Now demanding that this Hamiltonian involve the poten- 
tial V+\x) results in a differential equation that must be 
satisfied 

f( x ) + f( x )-Vi°\x)=0. 

This is a Riccati equation in its canonical form. The ex- 
plicit closed-form solution of this equation is not known 
typically, but one understands that the superpotential 
Wo(x) is a particular solution. This is enough to con- 
struct the general solution f(x) which depends on an 
arbitrary integration constant that can be considered as 
a free parameter in the partner Hamiltonian 



H = BiB = -di 



V?\x) 



2f'(x) = -dl + V(x). 



According to SUSY QM the potentials uj 0) (x) and V(x) 
are isospectral (except for the lowest state of V(x)) 
provided that f(x) is nonsingular. In addition, since 

BB^ = AqAq, it follows that the potentials V_ (x) and 
V(x) have strictly identical spectra. 

In ref. [ij Mielnik performed factorization of the har- 
monic oscillator potential in this manner. Mielnik ob- 
tained one-parameter family of potentials with the oscil- 
lator spectrum, but as we have just seen the procedure is 
straightforwardly generalized to any potential V_[°^ (x) . 

In the standard (i.e. based on the first-order intertwin- 
ing relation (1)) unbroken SUSY QM it is impossible to 
use an excited state of the original potential and at the 
same time avoid creating singularities in the partner po- 
tential @ . There is no guarantee that the resulting wave 
functions are normalizable and energy levels degenerate. 
The purpose of the present article is to modify the opera- 
tors B and B^ in such a way as to determine new strictly 
isospectral potentials without being forced to solve Ric- 
cati equations (by reducing the Riccati equation whose 
appearance in the factorization problems is typical to the 
solvable Bernoulli equation) and, more importantly, by 
applying the non-uniqueness of factorization to the su- 
perpotentials generated by the excited states of a poten- 
tial, since these also satisfy the Schrodinger equation. 



II. MODIFIED FACTORIZATION 

In this section we show the consequences of the non- 
uniqueness of factorization method extended to the ex- 
cited states of a potential, rather than just the ground 
state. In the literature the Hamiltonians and 



are called "bosonic" and "fermionic" respectively. We 
show that the degeneracy of energy levels of partner po- 
tentials depends on whether the bosonic or fermionic 
Hamiltonians admit non-unique factorization. 



A. Bosonic Hamiltonian 

Let there be given an analytically solvable non-singular 
potential (x) whose energy eigenvalues E n and wave 
functions ip n (x) are known. Without loss of generality, 

let E be zero, so that vd 0) (x) = ip%(x)/ip (x) = W£(x)~- 
Wq(x) and also define 

V M (x) = M{x)/M*) = - W^x) 

where W n {x) — —d x lmp n (x) is taken to be the super- 
potential corresponding to ip n {x). From the Schrodinger 



equation it follows that V^ n \x) — Vl v, (x) — E n , so that 

the potentials V_ (x) are non-singular, even though the 
superpotentials W n (x) are always singular for n > 0. Ad- 
justing the energy scale seems appropriate: one simply 
subtracts from the potential the energy of the excited 
state so that the resulting potential can be factored. 
Next we introduce the operators 

B n =d x + f(x) + W n (x) and B n = -d x + f(x) + W n (x) 

where f{x) will be determined below. Notice when n = 
these definitions reduce to the familiar case of standard 
unbroken SUSY QM if f(x) = and to the Mielnik's 
factorization Q if f(x) ^= 0. 

The factorization of the Hamiltonian 
leads to 



(0), 



BtB n 



H 



in) 



d 2 x + V M {x) + f{x) + 2W n {x)f(x) - f'(x). 



If we require that f 2 (x) + 2W n (x)f(x) - f(x) = 
the Hamiltonian becomes trivial because the potential 
U W (x) is related to (x) by a constant shift. On the 

other hand, the partner Hamiltonian = BnB^ is 

less trivial 



H 



(n) 



-dl + V^+2f{x) 



where uj n) ' (x) = W%(x) + W^x). The function f(x) is 
not arbitrary - it is a solution of the Bernoulli equation 
(a specific example of the Riccati equation) : 



f(x)=f(x)+ 2W n (x)f(x) 



and reads 



fn(x) = 



C-! X x Jn 2 {s)d S 



where C, xq are constants. It follows that ^ n {x) must be 
inverse square integrable; however, in general the wave 
functions do not possess this property. 
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There is yet another problem, namely, singularity of 
the potentials V_[ n ' (a;) for n^O corresponding to the ze- 
ros of the wave functions. Consequently, the breakdown 
of the degeneracy of energy levels of the Hamiltonians 

i n) and H ( + ] occurs (in addition to H_ and H^). 



B. Fermionic Hamiltonian 

The difficulties of establishing the degeneracy theorem 
for bosonic Hamiltonians suggest to reverse the order of 
the operators B n and B\ and start with the fermionic 
Hamiltonian HY 1 — 



H 



(n) 



dl + Vi n \x) + f(x) + 2W n (x)f(x) + f(x) 



where V±^ (x) are dehned as usual. We again obtain the 
Bernoulli equation 



f(x)+f(x)+2f(x)W n (x) 
whose general solution is 







fn{x) 



(2) 



where C, xq are constants and ip n (x) is assumed to be 
square-integrable. 

If it is possible to restrict the domain of the parameter 
C and make f n (x) free of singularities, then the potential 
V M (x) in 



= Bl B n = -dl 



y{n) = _Q2 



V 



(n) 



constitute a one-parameter family of potentials isospec- 
tral to the potential vi™** (x) . 

To see this note that the Schrodinger equation 
Hr'ipk = (E k -E n )ip k implies 

H^[BlA n ip k ] = B\B n B\A n i, k 
= BlA n AlA n ip k 
= (E k - E n )[B\A n ^ k ] 

where we have used the non-uniqueness of factorization 
of the Hamiltonian j n) = A n A\ = B n B\ L . So if i/> k (x) 
is an eigenfunction of the Hamiltonian H_ with en- 
ergy eigenvalue E k — E n , then B^Anipk is an eigenfunc- 
tion of H_ with the same energy. Similarly, from the 
Schrodinger equation H^-tp^ = E^-tp^ (where in 
Efc , k denotes the energy level and (n) refers to the n th 
eigenfunction of the Hamiltonian H_ ) it follows that 



H^\AiB n ^]=E^[AiB n ^} 



(n)r 



Hence, the normalized eigenfunctions of the Hamilto- 



nians H^' and are related by 



# n) (*0 



[E k - E n )- l [B\A n ^ k {x)\ 



(3) 



and 



E n y\AiB n ^\x)] 



where k 7^ n. The operators A n or B n destroy a node 
in the eigenfunctions, but they are followed respectively 
by the operators B^ or A\ that create an extra node. 
Thus, the overall number of the nodes does not change. 
In addition, the normalization does not require positive 
semi-definiteness of the energy eigenvalues, as in the stan- 
dard case. This is good because negative energy states 
appear when n > 0. 

For any n there is always one missing state k = n 
which can be obtained by solving the hrst order differen- 
tial equation = (by construction the state ij)^ 
has to be annihilated by the operator B n ): 



#(x) 



dx 



Therefore, 



= - W n (x) + 



c+rm S )ds 



ch- 



in 



ipr, 



c + r x ji{s)d S 



^(x)=N(C)x 



1pn 



~ (n) 

with the corresponding energy En — 0. All other energy 



(4) 



eigenvalues satisfy E 



(n) _ 



E k — E n . The normalization 



constant N(C) depends on the parameter C and other 
parameters of the potential such as width, depth etc. It 
is a constraint that allows one to determine the values of 
C for which the potentials V_ (x) are nonsingular and 

eigenfunctions ^^(x) are well-defined. 

One observes that the intertwining relationship be- 
tween the Hamiltonians and is of the second 
order: 

H ( l l) BiA n = BlA n H { - ] and A* n B n = 4 B n H™ . 

In the second-order SUSY QM 7} two different Hamil- 
tonians are intertwined by an operator of the second- 
order in derivatives, say, A = d% + r)(x)d x + 7(2;). If A 
can be written as a product of two first-order differential 
operators with real superpotentials, then we call it re- 
ducible (otherwise one refers to it as irreducible). Thus, 
our construction is equivalent to the second-order SUSY 
QM with the reducible operator A = — B^An. Perform- 
ing an explicit factorization one finds that —i](x) = f n (x) 

and —7(2:) = V_(x) + f n (x)W n (x). Pros and cons of 
these related approaches are discussed in detail in the 
concluding section. 
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From now on we will discuss the degeneracy of energy 

(n) ~ (n) 

levels of the Hamiltonians H_ and H_ only, leaving 

aside the Hamiltonian H_l' which plays an intermediate 
role in this construction. 



III. EXAMPLES 

Here we illustrate the results developed in the preced- 
ing section by providing examples that arise from well- 
known potentials and obtain some previously unreported 
potentials which might be of interest in various fields of 
physics and chemistry. One can also consult the ref. Q 
where factorizations of the harmonic oscillator potential 
were performed. 



A. Morse potential 

Let us first consider the Morse potential 



hi 



vH >) (x)=A i -B(2A + a)e 



i'X _j_ J^l ^ — IolX 



where the constants A, B and a are nonnegative. There is 
a finite number of energy levels Ek — ka{2A — ka) where 
k takes integer values from zero to the greatest value for 
which ka < A. For concreteness let us take A = 2 and 
a = B = 1. The partner potential V_ (x) is obtained 

from the ground state wave function ipo(%) = e~ 2x ~ e 
of the potential vl 0) (x) = 4 - 5e~ x + e~ 2x ; 



V { °\x) =4-5e~* + e 



-2x 



-4x-2e 



- 16 



dx \ C + e- 2e_x (3 + 6e~ x + 6e~ 2x + 4e~ 3:E ) 



As the potential V_ (x) it has only two bound states 

with eigenvalues Eq ^ = and E^ — 3. The normalized 
ground state wave function is 



8C(C+3) -2x-e~ 
3 e 



C + e- 2e ~* (3 + 6e- x + 6e~ 2x + Ae~ 3x ) 

Hence, the potential V_ (x) is nonsingular as long as 
C£ [-3,0] (see Fig. 1). 

The normalized wave function ^ (a;) is determined 
by applying the operator B^Aq to the first (and only) 
normalized excited state -01 (x) — 2/\ /r 3e~ x ~ e (3— 2e~ x ) 
of the potential V^\x): 

7(0) , . 2e- e ~*(6+l2e x + 9e 2x ) + Ce e ~* : (3e 2x -2e x ) 
Vi (x) =- 



\/3(4 + 6e x + 6e 2x + 3e 3x + C< 



We would like to remind the ladder operators for 
the wave functions of the Morse potential given in (5) 




C = -3.5 



FIG. 1: A few members of the one-parameter family of poten- 



(5) tials V_(x) isospectral to the Morse potential V}_ u> (x) with 



(0) , 



A — 2 and a = B = 1 (thick blue line) 

and explicitly derive them for the wave functions of the 
isospectral partner potential. Let's denote s = A/a and 
y = 2B/ae~ ax which is the common choice in the SUSY 
QM literature. Then for the creation K + and annihila- 
tion K_ operators we have Q: 



and 



K + = 



1/2 



2s - 2n - 1 



d u 



1/2 



2s- 2n + 1 



(we note that K- ^ KV) with the following effect 
K + ipk(y) oc tpk+i(y) and K-ipk+i(y) oc ipk{y)- The pro- 
portionality factors can be calculated after normalizing 



the eigenfunctions ipk(y) 



~ y / 2 L 2 f~ 2k (y) where 
L 2 , s ~ 2k (y) are associated Laguerre polynomials. 

The equation (3) enables us to deduce the ladder 

operators for the eigenvectors ipk(y) of the potential 
V_(x) whose energy spectrum is identical to that of 
the Morse potential V_ (x). The corresponding rais- 
ing and lowering operators for ipY 1 ' (y) with k ^ n are 
(BlA n )K+(AlB n ) and (B^A n )K- (A^B n ). Exploration 
of the higher-order ladder operators is the direct conse- 
quence of extending the first-order SUSY QM. 



B. CPRS potential 

In ref. ® Carinena, Perelomov, Ranada and 
Santander (CPRS) have studied the following one- 



5 



dimensional non-polynomial exactly solvable potential 
(we define our Hamiltonian to be — 2i?cRPS + 3): 

2x 2 



V W (x) = x 2 + 3 + , 



1 



(2x 2 



1) 



This potential asymptotically behaves like a simple har- 
monic oscillator but its minimum at the origin is much 
deeper than in case of the harmonic oscillator. Using 
SUSY QM techniques it was shown by Fellows and Smith 
10] that V^_°\x) is a partner potential of the harmonic 
oscillator x 2 + 5 and, therefore, their energy levels are the 
same. Here we further analyze the CPRS potential and 
find new potentials with the oscillator spectrum (see also 
ref. 1). 

The ground state energy Eq = and wave function 



e- x I 2 
2x 2 + 1 



of the potential (x) allows one to find its isospectral 
partner 



Vl v '(x) =x 2 + 'i + . 



(2x 2 + l) 2 



-8 A 



dx \ 2x(2x 2 + l)e- x2 + (2x 2 + 1) 2 (C + 0Ferf x) I 

which has no singularities when |C| > y/n (see Fig. 2) as 
follows from normalizing the ground state wave function 
-0q°^(x) (see below). 

Its eigenvalues are the same as that of the potential 



Ui 0) (x) and given by E%> = 2k + A for k 
normalized ground state wave function 



1,2 The 



^ 0) (x) - 



^2(C 2 -n)/^e-* 2 / 2 
2xe x2 + (2x 2 + 1)(C + y/netf(x)) 



corresponds to the energy eigenvalue Eq ^ = 0. The rest 
of the cigcnfunctions can be derived using equation (3). 

Neither Carinena et al., nor Fellows and Smith pro- 
vided the raising and lowering operators for the wave 
functions ipk (x) of the CPRS potential. Here we address 
the question of finding ladder operators for the CPRS po- 
tential and its isospectral partner. Taking into account 
that the CPRS potential itself is a partner of the har- 
monic oscillator, we obtain its raising and lower- 
ing A^aA operators where 



A = d x + x 



Ax 



2x 2 



1 



is needed to move between the CPRS potential and har- 
monic oscillator whose creation and annihilation opera- 
tors are a* and a respectively. Thus, the ladder opera- 
tors for the wave functions (x) of the potential ui™' 
become (Bl A n )A^ A(A\ l B Tl 
for k =/= n. 
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FIG. 2: Plot of the potential Ui 0) (x) with C = 1.8 (close to 
^/tt) and the unnormalized probability densities (dashed line 
at the corresponding level position) for its three lowest energy 
levels. The limit C — > oo corresponds to the CPRS potential 
(thick blue line). 



C. Infinite square well potential 

Despite its simplicity, the one-dimensional infinite 
square well potential with a deformed bottom requires 
some new techniques for obtaining solutions of the corre- 
sponding Schrodinger equation and usually one is unable 
to solve it exactly. In a recent paper [Tljj , exact solu- 
tion for the problem with sinusoidal bottom has been 
deduced. In this subsection we explicitly find potentials 
with undulating bottom and energy spectrum coinciding 
with that of the infinite square well. 

The wave functions and energy eigenfunctions of the 
infinite square well potential V_ (x) = —n 2 /L 2 of width 
L are given by ipk(x) = sin((fc + l)7ra;/L) with < x < L 
and Ek = k(k + 2)ir 2 /L 2 . Using this time for diverseness 
the first excited state wave function tpi (x) we find a pair 
of partner potentials, namely, the infinite square well po- 
tential with flat bottom 

V W (x) = -An 2 /L 2 

and the infinite square well potential with non-flat bot- 
tom also defined in the region < x < L (see Fig. 3): 



47T 2 

L 2 



16 d 

dx 



(2nx/L) 



C + Ax - L/ir sin (Anx/L) 



Both of the potentials have identical energy spectra 
j and {BlA n )A^aA{AlB n ) ^ = (fc - l)(fc + 3)tt 2 /L 2 . The normalized first ex- 
cited state of the potential V^}\x) is calculated from (4) 
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-10 - 



FIG. 3: Selected members of the family of one-parameter 
potentials V_(x). The limit C — > oo corresponds to the 



infinite square well V[ 0) {x) 
line) . 



and reads 



-4 of width L = 7r (thick red 



sin (2irx/L) 



C + Ax - L/ir sin (Anx /L) 



provided that C £ [— 4L, 0]. The wave functions ip(^\ x )i 
^>^\x), . . . can be found from (3). We only calculate the 
normalized lowest state eigenfunction: 

$1) {x) = ^f (MC + 4x) - 8L sin ^ + L sin *p) 



-3ir 2 /L 2 



3 J f (I, sin 3fE - 7r(C + 4tt)) 



(i) 



It corresponds to the negative energy Eft 

as expected since the potential V_ (x) is generated by 
the first excited state of the original potential. Note that 
the potential V_ (x) satisfies 

VP (C, x) = (C + 2L,x + L/2). 

It is known |8| that the eigenvectors tpk{ x ) of the 
Hamiltonian admit the following creation and an- 

nihilation operators: 



Mi 



and 



M_ = 



' 7TX \ 

cos | — )k 



cos ( — 1 k sin f— ^ 9. 



L 



sin 



fc _1 (fc- 1) 



where one defines the "number" operator k and its 
inverse k~ l such that kip^(x) = (k + l)ijjk{x) and 
k~ 1 ipk(x) — (k + l)~ lr ipk(x)- The ladder operators M± 
obey 

M--ip k (x) = kipk-i(x) and M + tp h (x) = (k + l)ip k +i(x). 
It is not hard to convince yourself that the raising 



of 



and lowering operators for the wave functions ip k 

~ (n) 

the partner isospectral Hamiltonian H_ are given by 
(BlA n )M+(AlB n ) and (fit A n )M-{Al B n ) respectively 
for k ^ n (when k = n use equation (4)). 



D. Two-parameter set of potentials isospectral to 
the harmonic oscillator 

Given an eigenfunction tp n (x) of the potential V^\x) 
one can find the wave function i/^™" 1 (x) of the one- 
parameter potential using the equation (3). Now 
one can repeat this procedure and instead of the eigen- 
function tp n (x) in (2) and (3) use ^ l \x) to obtain a 
two-parameter potential V^" ,fe ) (x) and its eigenfunctions. 
One can go on with this construction and obtain well 
defined multi-parameter potentials strictly isospectral to 
the potential v}_ k \x). 



-1 



Let us focus on the harmonic oscillator (x) = 
(with uj = 2) whose ground state wave functions is 

tpo(x) — e~ x I 2 . The potential V_ (x) is carefully dis- 
cussed in refs. 0, 0, EH each using different approaches, 
so in the following we omit unnecessary calculations and 
only state its normalized first excited state wave function: 



i>[ 0) (x) = 




3a;2 / 2 (l + 2Cxe x2 + ^xe x2 erf (a)) 



2C 



/nerf(x) 



where \C\ > \pn : /2 to guarantee non-singularity of the 
potential V^\x). Applying (2) to the wave function 
^^(x) we get the two-parameter potential (see Fig. 4): 



V^\x) 



dx\C + 0?/2erf(x) q + (^ 0) (s)) 2 ds J 

which is isospectral to the potential V_\x) — 2, which is 
in turn isospectral to the harmonic oscillator V^\x) — 
x 2 — 3, i.e. its energy levels are Ef' 1 ^ = 2(k — 1). 

The potential V_ (x) is non-singular for any C ^ 
and | C + 1/(4C)| > \/tt/A as follows from normalizing 
its ground state wave function. This family includes the 
oscillator potential x 2 — 3 in the limit C, C — > oo; the po- 
tential V^\x) arises when C — > oo; and finally V^'^(x) 
reduces to the potential V_ (x) Q in the limit C — >• oo. 
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FIG. 4: Plot of the potentials V"i 0,1) (a;), vP(x) = x 2 - 3 
and the non-normalized probability densities (dashed line at 
the corresponding level position) for the three lowest energy 
levels of V^' 1} (x). 



Let's briefly mention how to obtain its eigenfunctions. 
There is an expression similar to (3) for k = 2, 3, ... : 

^ 1] {x) cx B\Ai$p{x) cx B\AiBlA il> k (x) 

where (x) and ^ (x) are the eigenfunctions of the po- 
tential V_(x) and the harmonic oscillator accordingly. 
The operators B[, A\ are defined by 



A 1 =d x - d x In V^ 0) (x) 



and 



d x + d x In 



^i 0) (^) 



Lastly, the raising and lowering operators for the eigen- 
vectors "i/^'O) are given by S]JliBJA AJaJb ^i-Bi 
and AiB^oAoAjBo^i-Bi with Aj, A being the cre- 
ation and annihilation operators of the harmonic oscilla- 
tor. 

The two-parameter family of potentials with oscillator 
spectrum was also derived by the so-called second order 



intertwining technique in [12[. The advantages of the 
presented technique of getting multi-parameter sets of 
isospectral potentials are apparent. 



IV. CONCLUSION 

After the discovery of supersymmetry in string the- 
ory and then field theory, factorization was recognized 
as the application of supersymmetry to quantum me- 
chanics. The non- uniqueness of factorization serves as 
an avenue for the construction of many isospectral po- 
tentials. In this paper, we have explored the generality 
of this method by extending it to the excited states of a 
potential. Some nonsingular isospectral potentials that 
arise from the technique have been presented in this pa- 
per. These include one-parameter extensions of the well- 
known infinite square well and Morse potentials as well as 
not so familiar CPRS potential and two-parameter exten- 
sion of the harmonic oscillator. For some potentials the 
associated wave functions and probability densities have 
been derived and plotted. The ladder operators were 
determined explicitly. The application of this technique 
may be of significant interest because it can be applied 
to any one-dimensional quantum mechanical potential. 

The most general approach in the second-order SUSY 
QM is based on an arbitrary solution of the Schrodinger 
equation for the initial potential, rather than on its 
ground or excited state wave functions as discussed in the 
present article. However, there are certain advantages 
in such a presentation. For example, one can explic- 
itly construct the ladder operators for both isospectral 
Hamiltonians. It is also possible to avoid some technical 
complexities of the most general approach by mimicking 
the traditional first-order SUSY QM. For instance, in the 
second-order SUSY QM none of the expressions AA' or 
A^A coincide with any of the isospectral partner Hamilto- 
nians, but are quadratic forms in them. For comparison 
in our construction, which is based on the non-uniqueness 
of factorization, the appearance of the atypical Hamilto- 

(n) 

nian if_j_ at the intermediate stage does not affect the 



isospectral partner Hamiltonians H_ and 



r(n) 
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